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C\J ■ ABSTRACT 



Aims. Photometric data of galaxies covering the rest-frame wavelength range from far-UV to far-IR make it possible to derive galaxy 
■ properties with a high reliability by fitting the attenuated stellar emission and the related dust emission at the same time. 

Methods. For this purpose we wrote the code CIGALE (Code Investigating GALaxy Emission) that uses model spectra composed 
of the Maraston (or PEGASE) stellar population models, synthetic attenuation functions based on a modified Calzetti law, spectral 
(~ l , line templates, the Dale & Helou dust emission models, and optional spectral templates of obscured AGN. Depending on the input 

^ I' redshifts, filter fluxes are computed for the model set and compared to the galaxy photometry by carrying out a Bayesian-like analysis. 

CIGALE was tested by analysing 39 nearby galaxies selected from SINGS. The reliability of the different model parameters was 
evaluated by studying the resulting expectation values and their standard deviations in relation to the input model grid. Moreover, the 
influence of the filter set and the quality of photometric data on the code results was estimated. 



o 

^ . Results. For up to 17 filters with effective wavelengths between 0.15 and 160 yum, we find robust results for the mass, star formation 

' rate, effective age of the stellar population at 4000 A, bolometric luminosity, luminosity absorbed by dust, and attenuation in the 

far-UV. Details of the star formation history (excepting the burst fraction) and the shape of the attenuation curve are difficult to 
investigate with the available broad-band UV and optical photometric data. A study of the mutual relations between the reliable 
^ ' properties confirms the dependence of star formation activity on morphology in the local Universe and indicates a significant drop in 

' this activity at about 10" M© towards higher total stellar masses. The dustiest galaxies in the SINGS sample are present in the same 

I mass range. On the other hand, the far-UV attenuation of our sample galaxies does not appear to show a significant dependence on 

, star formation activity. 

• Conclusions. The results for our SINGS test sample demonstrate that CIGALE can be a valuable tool for studying basic properties of 

• ' galaxies in the near and distant Universe if UV-to-IR data are available. 

G^ ■ 
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• '"j 1. Introduction ticular, studies of broad-band SEDs suffered from difficulties to 

, . • , , disentangle age, metallic! ty, and dust effects. The lack of infor- 

■ Apart from dark matter, galaxies consist of stars, gas, and dust. ^^^-^^ ^^^^ necessary to assume a typical shape of the at- 

. Stars produce the radiation that allows us to observe galaxies ^^^^^^^^^ (Calzetti et al. 1994, 2000) and to apply simple 

in a wide wavelength range. Interstellar gas mainly modifies or optical to estimate the amount of attenua- 

the spectral energy distributions (SEDs) of galaxies by addi- ^-^^ p^^^^^j^ variations in the galaxy dust properties could not 

tional hne einission or absorption. Finally, interstellar dust af- considered in this way, which caused high uncertainties in 

fects galaxy SEDs by extinction, i.e. absorption and scattering .^i^eters such as the SFR. However, in view of many high- 

of stellar radiation, especially in the UV and optical and re- j^. ^^^^ collected by Spitzer, it is now possible to study 

emission of the absorbed energy in the IR The dust-induced ^^^^ j ^ j^. ^.^^^j^.^^^ ^ ^.^^ 

energy conversion makes it difficult to study basic properties re- ^^^glength range, which allows us to better understand dust ef- 
lated to the stellar populations like the star formation rate (SFR), f^^i^ 
since a good knowledge of the galaxy SEDs over a wide wave- 
length range is necessary to evaluate the effect of dust. In face of the complex SEDs of galaxies consisting of emis- 
In the past, the availability of IR galaxy data (IRAS and ISO) sion from different stellar populations which is modified by in- 
was mainly restricted to the nearby Universe. Hence, most low- terstellar gas and dust in a complicated way, it is clear that the 
to-high-redshift galaxies could be observed in the rest-frame UV derivation of galaxy properties needs a reahstic modelHng of 
to near-IR only. Since the total amount of dust emission in the galaxy SEDs. For the stellar populations those models were pro- 
IR could not be used to determine the amount of attenuation of duced by, e.g., Fioc & Rocca-Volmerange 019970 . Bruzual & 
the stellar continuum in the observed wavelength ranges, in par- Chariot (2003), and Maraston (2005 ). Dust attenuation curves 

were either derived by studying the SEDs of nearby star-forming 

Send offprint requests to: S. Noll galaxies (Calzetti et al. 1 1994] [20001 ) or by more theoretical radia- 
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live transfer models (e.g., Witt & Gordon |2000| ). The emission 
of dust grains in the IR and polycycUc aromatic hydrocarbons 
(PAHs; see Peeters et al. 2004 and references therein) in the mid- 
IR was described by templates of Chary & Elbaz ( 120011 ). Dale & 
Helou (|2002]| Lagache et al. (120031 120041) . and Siebenmorgen & 
Rriigel (120071) . Models that also include stellar emission and dust 
attenuation at shorter wavelengths were published by, e.g., Silva 
et al. ( [1998] ), Dopita et al. (2005 ), and da Cunha et al. (2008 ). The 
former also consider the evolution of dust properties dependent 
on the age of the stellar population. Sophisticated fitting codes 
of galaxy SEDs mainly focus on the derivation of redshifts and 
relatively robust galaxy properties like masses (e.g., Bolzonella 
et al. |20D0; Feldmann et al. 2006; Walcher et al. .2008j . Models 
for the IR are usually not included. 

The above collection of publications on the modelling of 
galaxy SEDs shows that the comparison of data and models is an 
important tool for studying the physical properties of galaxies. 
However, details of the listed models like galaxy types, wave- 
length regimes, and characteristic parameters can differ a lot. 
On the other hand, the different astronomical questions and the 
available data sets require suitable models. Finally, optimised 
routines are necessary in order to derive the characteristic galaxy 
properties under investigation from the comparison of data and 
models. 

Hence, we feel that an observer-friendly fitting code for star- 
forming galaxies at different (given) redshifts is still missing 
which especially provides star formation histories, dust attenua- 
tion properties, and masses and uses photometric data covering 
the rest-frame UV to IR to constrain these properties. In other 
words, we are interested in a code that calculates the effect of 
dust on galaxy SEDs in a consistent way in order to obtain prop- 
erties (such as the SFR or the effective age of the stellar pop- 
ulations) which are important for understanding the evolution 
of galaxies. Based on a procedure described in Burgarella et al. 
(120051 ). we have developed a code that derives galaxy properties 
by means of a Bayesian-like analysis. In Sect. |2] we describe the 
models used and the fitting procedure. In Sect.[3]we test the code 
for nearby galaxies selected from the SINGS sample (Kennicutt 
et al. 2003 ) for which good photometric data from the UV to the 
IR are available. Finally, the resulting findings concerning the 
applicability and effectivity of the code are discussed in Sect. |4] 



2. The code 

Our programme package CIGALI^ (Code /nvestigating GALaxy 
Emission) is characterised by a series of modules which are fed 
by spectra and parameter files. A flow chart of the code is shown 
in Fig.[T] The general idea is to build stellar population models 
first. Secondly, the dust is considered by reddening the stellar 
SEDs and re-emitting the absorbed energy in the IR. Dust emis- 
sion in the IR due to non-thermal sources can also be added. 
Interstellar lines are taken into account as corrections of the 
dust-affected SEDs. Redshift-dependent filter fluxes of the en- 
tire model set are, then, calculated for a direct comparison to the 
input object data. Finally, depending on the for the best-fit 
models of a set of bins in the parameter space, probability distri- 
butions as a function of the parameter value are calculated and 
used to derive weighted galaxy properties. Before we discuss 
this Bayesian-like analysis and the interpretation of its results in 
more detail in Sect. 12.21 and 12.31 we start with a description of 
the models used in Sect. 12.11 



' Our code is provided at http: //www. oamp . fr/cigale/. 



2.1. Models 

For CIGALE we take models indicating the net emission from 
a galaxy in the wavelength range from far-UV to far-IR. This 
means stellar emission, absorption and emission by dust, and at 
least the strongest interstellar absorption and emission lines have 
to be considered. 

2.1.1. Stars 

Concerning the stellar SEDs we have decided to focus on the 
models of Maraston (2005), since they consider the thermally 
pulsating asymptotic giant branch (TP-AGB) starfl in a real- 
istic way. These stars are bright intermediate-age stars of 0.2 
to 1-2 Gyr which mostly contribute from red optical to near- 
IR wavelengths. Therefore, they are particular important for a 
reliable stellar mass determination, but star-formation-related 
parameters are also affected. The insuflicient consideration of 
TP-AGB stars in older but widespread models of Bruzual & 
Chariot (2003T) and Fioc & Rocca-Volmerange (fT997 ; PEGASE) 
typically increases the mass by 0.2 dex for star-forming galax- 
ies with an important population of intermediate-age star^ 
(Maraston et al. 2QQ6, Salimbeni et al. ,2009j . Therefore, 
CIGALE also allows PEGASE models to be used for backwards 
compatibility. 

For the calculation of complex stellar populations (CSPs) we 
consider single stellar populations (SSPs) of Maraston (2005) 
and PEGASE of different metallicity Z and Salpeter tl955i ) or 
Kroupa (2001) initial mass function (IMF). Masses based on the 
Salpeter IMF are about 0.2 dex higher than for the Kroupa IMF. 
The galaxy mass is normalised to 1 M© and comprises the total 
mass of the stars (active and dead) plus gas that originates from 
stellar mass loss. This means that all SSP models start with a 
stellar mass fraction of 100%. For a Salpeter IMF the fraction 
reaches about 70% at the age of the Universe. The total stel- 
lar mass is provided by the code. For obtaining CSPs the SSPs 
of different ages are weighted and added according to the star 
formation scenario chosen. We provide "box models" with con- 
stant star formation over a limited period and "t models" with 
exponentially decreasing SFR, respectively. While for the for- 
mer scenario and ongoing star formation the instantaneous SFR 
at look-back time f' = is simply calculated by the galaxy mass 
divided by the age, i.e. Mgai/f, the SFR of the latter scenario re- 
sults in 



where r is the e-folding or decay time. We focus on these r mod- 
els for the rest of the paper. 

CIGALE allows two CSP models to be combined. Both 
model components are, then, linked by their mass fraction. 
This scenario makes it possible to consider bursts on top of an 
older passive or more quiescent star-forming stellar population. 
Although this approach allows galaxy SEDs to be fit quite well, 
it is clear that this scenario can only be a rough representation of 
the real star formation history (SFH) of a galaxy. On the other 
hand, the more complex SFH scenarios are, the more SED fitting 



- The Maraston models include TP-AGB spectra up to 2.5 fim only. 
Therefore, the flux level in a narrow wavelength range in the near-IR 
can still be systematically underestimated. 

^ For early-type galaxies in the nearby Universe the difference is ex- 
pected to be lower, since most stars are distinctly older than the maxi- 
mum age of a TP-AGB star. 
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Fig. 1. Flow chart of CIGALE. The different programme modules and the corresponding parameter and template inputs are shown 
from the left to the right. The alternative (and less preferred) use of the PEGASE models is marked by a dashed arrow. The other 
dashed arrows refer to the option to take either IR photometry directly (consideration of Dale & Helou and possibly Siebenmorgen 
et al. models) or externally estimated IR luminosities for considering the IR regime. The presented flow chart assumes that the 
object redshifts are taken from the photometric input catalogue. 



suffers from degeneracies. In view of the uncertainties in the ba- 
sic stellar population model parameters, CIGALE also computes 
two effective parameters: the mass-weighted age 

r'f'SFR(f')df' 



and the age fD4000 derived from the D4000 brealfl (see Balogh 
et al. 11999) of the unreddened SED for a single stellar popu- 
lation, which shows a quasi-monotonic increase of D4000 with 
age (see Fig. |2]l. The two ages are complementary, since they 
trace the old mass-dominating and the young light-dominating 
stellar population, respectively. 

2.1.2. Dust attenuation 

The absorption and scattering of star light by dust is described by 
extinction curves. In the case of galaxies these laws are preferen- 
tially called attenuation curves because of the dust-induced re- 
scattering of radiation in direction towards the observer, which is 
negligible if only a single star without a circumstellar dust shell 
is taken into account (see Krugel |2009l l. For typical sightlines 
towards Milky Way stars the extinction curve is characterised 
by increasing extinction at shorter wavelengths, which causes 
a reddening of the spectra, and the so-called UV bump at about 
2175 A, which produces a broad absorption feature (e.g., Stecher 
[19691 Cardelh et al. [19891 Fitzpatrick & Massa [20D7] ). For typ- 
ical sightlines towards the Large and Small Magellanic Cloud 
(LMC and SMC) the extinction curves are steeper and the UV 
bumps are weaker than in the Milky Way (Gordon et al. 120031 ). 



ratio of the average flux per frequency unit of the wavelength ranges 
4000-4100 A and 3850-3950 A 




1 1.5 2 2.5 

D4000 



Fig. 2. Relation between the age fD40{)0 and the strength of the 
D4000 break for Maraston (2005) single stellar population mod- 
els with Salpeter IMF and metallicities 0.5, 1, and 2Zq. Local, 
low amplitude extrema in the curves were smoothed out to avoid 
any ambiguities. 

The most extreme curve is found for the SMC, where the 2175 A 
feature is almost vanished completely. A non-existing UV bump 
and a moderate far-UV rise characterises the effective, average 
attenuation curve of nearby starburst galaxies found by Calzetti 
et al. a 9941 [2000] ). However, about 30% of the UV-luminous 
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Fig. 3. Illustration of synthetic attenuation laws differing in cr 
and 6. Apart from a pure Calzetti law (cr = and 6-0; solid 
curve), Calzetti laws modified by multiplication of a power law 
are shown. Lower 6 cause steeper slopes in the UV. In two cases 
the change of Ry by the power law is compensated by cr - -6 
(dashed curves). 

star-forming galaxies at 1 < z < 2.5 seem to exhibit a signifi- 
cant 2175 A feature which can be as strong as those found in the 
LMC (Noll & Pierini 2005; Noll et al. 2007 2009). Moreover, 
the typical width of the UV bump of these galaxies is only about 
60% of those measured for characteristic Milky Way and LMC 
sightlines. 

These observations show that a single attenuation curve as 
given by Calzetti et al. (120001) and often used in the literature is 
probably insufficient to describe star-forming galaxies, in par- 
ticular, at higher redshifts. Hence, we apply a more complex 
ansatz to describe the attenuation in our model galaxies. Since 
the Calzetti law seems to provide reasonable amounts of attenua- 
tion in starburst galaxies at least as an average for large samples, 
we take this frequently used curve as basis. For wavelengths be- 
low 1200 A, for which the Calzetti law is not defined, we lin- 
early extrapolate by using the value and slope at 1200 A. This ap- 
proach, which resembles the method of Bolzonella et al. (.2000J , 
avoids extreme slopes at low wavelengths that are not supported 
by observations (Leitherer et al. ' 20021 ). 

At first, the Calzetti law kiA)+Rv,o (Rv,o - 4.05+0.8) is mod- 
ified by adding a UV bump which is modelled by a Lorentzian- 
like "Drude" profile 



2, ,2 



D 



(i) = 



(^2 _ ^„2-)2 + ^2y2 ' 



(3) 



where Aq, y, and /Sbump are the central wavelength, width 
(^ FWHM), and amplitude (or maximum height), respectively 
(Fitzpatrick & Massa , 19901 12007> Noll et al. l20091 l. In order to 
produce different slopes, the modified Calzetti law is multiplied 
by a power law (A/Av)^, where Ay - 5500 A is the reference 
wavelength of the V filter. This choice allows the slope to be 
changed without altering the scaling parameter visual attenua- 
tion Av- Slopes 6 > Q represent flatter slopes than the Calzetti 



law, while 5-0 exactly reproduces this basic curve if no bump 
is considered (see Fig. [3]). If the resulting curve gives a negative 
attenuation correction, which even happens for a pure Calzetti 
law in the near-IR, no correction is applied. Since 6 causes 
a change in the original scaling parameter of the Calzetti attenu- 
ation curve Eb-v, the real Ry = Av/Eb-v differs from the value 
Rv,Q = 4.05 of the Calzetti law. In order to compensate for this 
effect and to have more freedom in the choice of this constant, 
the attenuation law is multiplied by the factor (1 - flRCR)//?v,o- 
For Or = 1.12 this term guarantees roughly constant Ry if 
Cr + 5 = const (see Fig. [3]l. For reasonable slope corrections 
-0.3 < <5 < 0.3 the real Ry changes by about +1% or less. 
In summary, the dust attenuation in magnitudes as function of 
wavelength is 



A(A) = Av 



1-L12^ 
<• " 4.05 



^5500 A/ 



The attenuation correction is applied to both possible stellar 
population components (see Sect. 12.1. Il l individually by allow- 
ing for different Ay. In practice, the visual attenuation of the 
young T model Ay ^sp and a reduction factor of the attenuation 
for the old t model /att are used. This simple procedure takes into 
account that stars of different age can be affected by different 
amounts of dust (Silva et al. il998i Pierini et al. ,2004i Panuzzo 
et al. 120071 Noll et al. l2007l l. 

In view of the complexity of the attenuation curves and the 
dependence of the total attenuation on the SFH due to the intro- 
duction of /att, we characterise the effective obscuration of the 
stellar radiation by two additional parameters that are derived 
from the final model SEDs. Apuv and Ay are defined as the ef- 
fective attenuation factors in magnitudes at 1500 + 100 A and 
5500+ 100 A, respectively. Both parameters are complementary. 
While Apuv probes the dust obscuration of the young stellar pop- 
ulation only. Ay also traces the attenuation of older and cooler 
stars. The latter parameter equals Av ^sp if only one r model is 
used to describe the SFH of a galaxy. 

2.1 .3. Dust emission 

Dust re-radiates energy predominantly absorbed in the UV and 
optical at relatively long, mid-IR to sub-mm wavelengths. In 
the far-IR dust emission is characterised by continuum radia- 
tion that originates from dust grains in thermal equilibrium. The 
wavelength distribution of the emission depends on grain tem- 
perature, i.e. longer wavelengths indicate lower temperatures. In 
the mid-IR the continuum radiation of star-forming galaxies is 
mainly caused by stochastically-heatecfl very small grains (see, 
e.g., da Cunha et al. 120081 and references therein). This emis- 
sion is superimposed by prominent broad emission features of 
polycyclic aromatic hydrocarbon (PAH) molecule^ (see Puget 
& Leger [19891 Sturm et al. l2000l Draine|2003l Peeters et al. 
12004 ). In contrast, mid-IR spectra dominated by hot dust contin- 
uum emission and lacking significant PAH emission are usually 
related to active galactic nuclei (AGN). 

The IR properties of star-forming galaxies can be described 
by physically motivated multi-parameter models (e.g., Dopita et 



' Thermal emission of very small grains in the mid-IR is restricted to 
extremely intense heating environments that are not typical of normal 
star-forming galaxies (e.g., Dale et al. 2001 1. 

* Very small grains and especially PAHs in their low energy states 
may significantly contribute to the emission at mm and cm wavelengths. 
In particular, spinning PAHs are proposed as the origin of the "anoma- 
lous emission" at these wavelengths (see Ysard & Verstraete 1 2009 i l . 
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Fig. 4. Illustration of a set of complete models differing in Ay 
(see legend). The higher Ay is, the lower the flux in the UV and 
the higher the flux in the IR is. All models are characterised by 
solar metallicity, t - lOGyr, f = 1 Gyr, Calzetti law (i.e. cr = 
and 6 = 0), an LMC-like UV bump, no AGN contamination, and 
a = 2.0. 



al. 120051 Siebenmorgen & Krugel|2002l da Cunha et al. l2008] l. 
However, since the IR data of galaxies at higher redshifts are 
usually quite sparse, we have decided to rely on semi-empirical, 
one-parameter models (Chary & Elbaz i2001i Dale & Helou 



a: 1.0 2.0 3.0 ~ 
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Fig. 5. Illustration of a set of complete models differing in a (see 
legend). The higher a is, the larger the peak wavelength of the 
dust emission is. All models are characterised by solar metallic- 
ity, T - 10 Gyr, t = 1 Gyr, Calzetti law (i.e. cr = and 6 - 0), an 
LMC-likeUVbump, no AGN contamination, and Ay - 1.2 mag. 



[20021 Lagache et al. |2003ll2004li. For CIGALE we take the 64 
templates of Dale & Helou ( 120021 ), which are parametrised by 
the power law slope of the dust mass distribution over heating 
intensity a. For higher a the contribution of relatively quiescent 
galactic regions characterised by weak radiation fields becomes 
more important and the dust emission peaks at longer wave- 
lengths (see Dale et al. 12001 ). The PAH emission pattern of the 
Dale & Helou (|2002]) models shows only little variation. This 
agrees with the observations as long as the IR emission is not 
significantly affected by an AGN (e.g.. Peelers et al. |2004| l. As 
for the Chary & Elbaz (2001) and Lagache et al. (120031 ) models, 
a characteristic IR luminosity Lir can be assigned to each Dale 
& Helou ( 120021 ) template. Those calibrations for luminosity- 
dependent IR SEDs of the Dale & Helou models are available 
from Chapman et al. (12003) and Marcillac et al. ( 120061 ). Since 
individual galaxies can significantly differ from these relations 
for sample averages, we do not directly use them for CIGALE. 
However, they can be taken to adapt the a parameter space for 
the kind of objects investigated. 

The semi-empirical Dale & Helou templates include stellar 
emission in the near-IR range and below (see Dale et al. 2001). 
In order to combine them with stellar population models the stel- 
lar contribution needs to be subtracted. Therefore, we scaled 
a 5Gyr-old passively-evolving Maraston model to the flux at 
2.5-3 jjm in the different Dale & Helou models to remove the 
stellar emission. For wavelengths below 2.5 //m we set the flux 
in the IR models to zero, in any case. The flux reduction drops 
below 50% at about 5 ^^m. The choice of the stellar population 
model is not critical, since the slope of the stellar continuum at 
these long wavelengths shows only little variation for different 
SFHs. 

The IR templates are linked to the attenuated stellar popu- 
lation models by the dust luminosity Ldust, i-e. the luminosity 
"absorbed" by the dust and re-emitted in the IR. Hence, the scal- 
ing of the Dale & Helou templates is not a free parameter. The 
change of the luminosity transfer with increasing Ay for a fixed 
attenuation law is illustrated in Fig. |4] On the other hand. Fig. |5] 
demonstrates for fixed Ljust the change in the shape of the IR 
SED if the only IR-specific parameter a is modified. 

The equality between the dust-absorbed and dust-emitted lu- 
minosity can be violated by dust emission caused by a non- 
thermal source. In particular, highly dust-enshrouded AGN rep- 
resent a problem, since they are difficult to identify in the UV and 
optical. Consequently, galaxies with such an AGN contribution 
look like normal star-forming galaxies in the far-UV-to-far-IR 
wavelength range excepting a warm/hot dust emission in the IR. 
Therefore, we allow for an additional IR dust emission compo- 
nent which is not balanced by dust absorption of stellar emission 
at shorter wavelengths. 

Siebenmorgen et al. (|2004al |2004bl ) provide almost 1500 
AGN models differing in the luminosity of the non-thermal 
source, the outer radius of a spherical dust cloud covering the 
AGN, and the amount of attenuation in the visual caused by the 
cloud. In principle, all models can be fed into CIGALE. Focusing 
on SEDs providing PAH-free mid-IR emission, the number of 
suitable models significantly decreases, however As reference 
model we take L = 10'^ L©, R = 125 pc, and Ay - 32 mag. 
Using this or similar models allows us to disentangle IR dust 
emission components caused by stellar and AGN radiation, re- 
spectively, since the AGN-related component peaks at signifi- 
cantly shorter wavelengths than the stellar one. Admittedly, the 
real SEDs of obscured AGN could significantly deviate from our 
preferred model. However, this is almost impossible to test with 
broad-band photometric data only. 
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2.1.4. Gas 

Strong spectral lines or blends can easily change filter-averaged 
fluxes by 10% or more. Therefore, spectral features need to be 
considered to avoid systematic deviations in the spectra. The cru- 
cial features are nebular emission lines of, e.g., HI, [Oil], and 
[O III] in the optical and interstellar absorption lines of different 
ions in the UV. The UV is also affected by stellar wind lines like 
CIV and nebular emission, i.e. essentially the strongly-varying 
Lya line. The narrow nebular emission lines in the IR are not ex- 
plicitly considered, since we assume that they already contribute 
to the low -resolution semi-empirical templates of Dale & Helou 
(I2OO2I 1. 

The spectral line correction is performed by using empirical 
line templates. We have taken the Kinney et al. ( 1996) starburst 
spectra to derive two optical emission line templates mainly dif- 
fering in the strength of the oxygen lines. For the UV we have de- 
rived two line templates from composites of high-redshift galax- 
ies of Noll et al. ( |2004| l, since their UV spectra have a better 



quality than those of Kinney et al. (11996 ). The UV line tem- 
plates include all lines between Lyman limit and 3000 A ex- 
cepting those which are already present in the stellar population 
models. The selection of the optical and UV templates depends 
on the flux ratio of the wavelength ranges 1420 - 1490 A and 
3400 - 3600 A in the attenuated stellar population models in or- 
der to reproduce the change of line strengths with UV continuum 
reddening in star-forming galaxies. The best agreement with the 
observed spectra is achieved if for ratios > 3.7 the templates 
with strong emission and weak absorption features are taken. 
For ratios < 1.0 and star-related Ldust/i'boi < 0.1 the spectral line 
templates are not considered at all in order to reproduce mainly 
passively-evolving galaxies. Before the selected line templates 
are added to the models, the optical line templates are scaled 
to the average flux in the wavelength range 3400 - 3600 A. On 
the other hand, the UV line templates are, first, multiplied to the 
dust-attenuated models cleaned from stellar lines by interpola- 
tion. The effect of the two sets of UV and optical line templates 
on model spectra is visible in Fig.|4] 

The described procedure for the consideration of interstel- 
lar lines in the model SEDs is, of course, relatively simple. 
However, it avoids systematic errors, although with significant 
uncertainties for individual objects with SEDs deviating from 
the composites analysed. A clear advantage of this phenomeno- 
logical approach is certainly the avoiding of new free parameters 
in the code. 



2.2. Fitting and parameter probability distributions 

The observational input data of CIGALE are photometric filter 
fluxes. Hence, the filter fluxes for the entire set of models are de- 
rived for the known redshifts of the objects given. The final mod- 
els taken also include a correction of the redshift-dependent ab- 
sorption of the intergalactic medium (IGM) shortwards of Lyo'. 
Here we take the more recent algorithm of Meiksin (2006 ) in- 
stead of the more frequently used corrections of Madau ( 1995i) . 

For sparse IR data it is possible to provide an externally 
computed IR luminosity to the code instead of using the IR fil- 
ters directly. For example, the Dale & Helou (,2002j models can 
be used in combination with the calibrations of Chapman et al. 
( 120031 ) and MarciUac et al. ( 120061 ) to derive typical Lir from 
MIPS 24 /im flux densities. In this case, Lir is converted into 
flux units for an artificial filter centred at 100/zm. This is, then, 
compared to the dust luminosity Ldust derived for each model 
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Fig. 6. PDF derivation for the SFR of a test galaxy. About 
7 X lO*" models are shown (small dots). The parameter space 
is divided into 50 bins, which are indicated by histograms. The 
best model of each bin is marked by a filled circle. The upper 
panel exhibits the;^^ distribution as a function of SFR. The lower 
panel shows the corresponding probabilities for the normalised 
X^, i.e. "A'besf Expectation values and standard deviations for 
the SFR can be derived from the probability distribution of the 
bin-specific best-fit models. 

converted to a filter flux as well. In order to save computing time 
no IR models need to be considered in this case. 

The comparison of model and noise-afffected object photom- 
etrjQ, i.e. of /mod,/ (per Mq) and /obs,/ for k filters, is carried out 
for each model by the minimisation of 



X' 



(A/gai /jTiod,/ /obs,/) 



(5) 



obs,/ 



with the galaxy mass Mgai (in Mq) as a free parameter The sta- 
tistical photometric errors are considered by o-gbs,/. The resulting 
minimum x' of each model can be compared to determine the 
model showing the best x' of the entire model grid for the ob- 
ject SED investigated. Moreover, the model-related allow us 
to perform a more sophisticated analysis to obtain galaxy prop- 
erties based on probability distribution functions (PDFs). 



' Since filter-averaged fluxes are compared, the calibration of pos- 
sibly considered Spitzer MIPS flux densities has to be changed to 
vfy = const. Consequently, the observed flux densities at 24, 70, and 
160 yum have to be multiplied by 1.040, 1.070, and 1.043, respectively 
(cf. MIPS data handbook), before they can be used in the code. 
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Fig. 7. PDF derivation for Ay ^sp. In this example, the discrete 
model Ai/ysp values increase in steps of 0.15 mag. For the de- 
scription of the plot see Fig.|6] 



Probabilities for individual models are often computed using 
the exponential term e^^'^^ (e g-, KaufFmann et al. 12003a! Salim 
et al. 120071 Walcher et al. .2008) . Then, PDFs for each param- 
eter can be derived by calculating the probabiUty sums of the 
models in given bins of the parameter space. However, the re- 
sults of this 'sum' method depend on the model density in the 
parameter space. In the case of a bad choice of the input model 
parameter values, this can cause an unintentional bias. Although 
the introduction of a lower threshold probability for the consid- 
eration of models can alleviate this effect, we adopt a different 
approach based on the best-fit models for given bins in the pa- 
rameter space and integrated probabilities, which we describe in 
the following. We refer to this approach as 'max' method. We 
compare the results of both methods in Sect. 13.2.21 

The likelihood of a model can be inferred from by inte- 
grating the corresponding probability density function from 
to infinity: 



n^o _^(n/2)-lg-A/2 

^^''"^^ L 2«/2r(n/2) 



(6) 



F denotes the Gamma function and n - k - I the filter-related 
degrees of freedom. Since the minimum x^, i-e. the fit quality, 
can vary a lot for an object sample, we prefer to use normalised 
= ~ Xhest (although not forced by the code) for com- 



^norm 



probability distributions in the model sets for the different ob- 
jects. This means that the best-fit model always has p - I- 

For calculating the object-related expectation values and un- 
certainties for the different model parameters, PDFs P{x) de- 
pending on the parameter value x are necessary. This requires 
to link the p{x^, n) for all m models studied to the P(x). We do 
so by introducing a fixed number b of equally-sized bins for each 
parameter. The range of bins is delimited by the lowest and high- 
est value of a parameter in the model set. We now derive the 
characteristic probability P,- of each bin / by searching the maxi- 
mum Pjix^, n) of the models j located inside the bin /. The prior 
Uji = 1 if a model j belongs to the bin /, otherwise ay, = 0. In 
summary, this procedure can be written in a Bayesian-like style: 



max ipjix ,n)aji). 

j=l,...,m 



(7) 



The resulting PDF Pi(x) envelops the distribution of models in 
the x-p plane (see Fig. |6]for a continuum of x and Fig. [7] for 
discrete x). The already mentioned advantage of this approach is 
that Pi(x) does not depend on the model density as a function of 
X. Consequently, the arbitrary choice of values for a parameter 
can hardly change Piix) as long as the covered range of x is 
comparable. 

Taking the Piix) as weights for each bin, the expectation 
value of each parameter is given as 



(x) 



(8) 



For Xj we directly use the parameter value of the best-fit model 
of each bin. It corresponds to the mean value of all models in 
a given bin if the number of bins is significantly higher than 
the number of realised parameter values, which should be the 
case for most model parameters. Finally, the standard deviation 
is derived by 



cr, = 



(9) 



CIGALE provides results for 16 basic input parameter^ the 
scaling parameter Mgai, and 10 additional output parameters that 
depend on different basic model properties. A complete list is 
given in Table [T] Apart from mean value and error, the code out- 
puts the parameter values of the best-fit model and the full PDFs 
for each parameter if desired. 

2.3. Interpretation of the code results 

The mean values (x) and errors cr^- derived by CIGALE have to 
be taken with care. The high complexity of the models and the 
limited number of photometric filters can cause that part of the 
parameters are not well constrained. Moreover, restrictions of 
computing time and hard-disc space limit and/or inherent limi- 
tations of the models used restrict the range and number of pos- 
sible values for each parameter which can be investigated by the 
code. Since good a-priori estimates of the PDF of a parameter are 
rare, the selected model grid may fail to cover the entire set of 
parameter values consistent with the data. In this case the calcu- 
lated average is probably misplaced and the errors are too small 



puting p in order to guarantee similar (and better comparable) 



The true number of simultaneously analysed, basic parameters is 
lower than 16 and should not exceed 10 in most cases, since parameters 
such as the metallicity or the central wavelength of the UV bump are 
usually defined by a single value. 



8 



S. Noll et al.: Analysis of galaxy SEDs from far-UV to far-IR with CIGALE: Studying a SINGS test sample 




Table 1. Description of the output parameters of CIGALE 



Fig. 8. Illustration of the change of the mean values and standard 
deviations of a Gaussian PDF by the limitation of the parameter 
range. The cut probability distribution (filled area) indicates a 
mean closer to the centre of the parameter range and an error 
smaller than the true value. 
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Fig. 9. Deviation of the measured from the real mean values and 
errors for a Gaussian PDF. The axes show the apparent expec- 
tation values and standard deviations scaled to the half param- 
eter range. The grid of solid curves indicates the true (and also 
scaled) mean values and errors for steps of 0.05 from to 1. 
Curves for constant real mean values -0.3, -1, and -3 outside 
the covered parameter range and variable real cr are shown by 
dashed lines. The solid curve crossing the grid delimits the trust- 
worthy area where the deviation between apparent and true mean 
values and errors is lower than 1% and 10%, respectively. 
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(see Fig.[8]l. Therefore, we study the uncertainty and reliability 
of the results of CIGALE by means of simulations for different 
PDFs for different sets of parameter values. 

At first, we consider a Gaussian PDF for an infinite number 
of values. By changing the position and width of the Gaussian in 
a fixed parameter range of normalised width 2, we obtain the re- 
lation between apparent and true mean values and errors as pre- 
sented in Fig. |9] For the area below the intersecting solid curve 
the regular grid in steps of 0.05 does not change significantly, 
i.e. A(ji:) < 0.01 and Acr^ < 0.1. For the centre of the Gaussian 
closer to the margin of the covered range and higher widths of 
the Gaussian the differences between true and apparent expec- 
tation values and standard deviations increase more and more 
until a curve that delimits the permitted area is reached. Close to 
this enveloping curve the corrections are enormous and highly 
uncertain because of the high density of curves. We can there- 
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Fig. 10. Curves enveloping the area allowed for data in standard- 
ised diagrams showing mean values and errors scaled to the half 
parameter range. The lowest solid curve corresponds to the en- 
veloping curve of the grid shown in Fig.|9l i.e., this curve is for a 
Gaussian PDF with the true mean inside the parameter range 
and an infinite number of parameter values. The dash-dotted 
curve shows the same for a box-shaped distribution instead of 
a Gaussian. The results for double Gaussians with 2 cr and 3 cr 
separation are marked by the dashed lines. The change in the 
data-covered area by the reduction of the parameter input values 
is indicated by the solid lines labelled by the number of values. 

fore divide the area of a standardised diagnostic diagram in three 
sub-areas which characterise reliable, uncertain/unreliable, and 
impossible results. 

For data points that are not too close to the enveloping curve, 
the diagram could, in principle, be used to recover the true mean 
values and errors. However, the real situation is more compli- 
cated, since PDFs do not need to be Gaussian (see, e.g., Figs.|6] 
and|7]). Fig. [To] shows enveloping curves (for the true mean lo- 
cated inside the parameter range) for different shapes of the 
probability distribution. For a box-shaped distribution the possi- 
ble errors are generally lower. On the other hand, double-peak 
profiles can cause relatively high apparent errors close to the 
centre of the parameter range. In this case, the true mean can 
be located in the other half of the parameter range, which is not 
possible for the single-peak PDFs discussed before. 

Finally, Fig. [10] illustrates that a finite number of parame- 
ter input values causes an increase of the errors. In the most 
extreme case of two values, only the margins of the parameter 
range are covered by data points, which produces significantly 
higher standard deviations than for a smooth PDF. In contrast to 
uncertainties in the shape of the probability distribution that im- 
pede a detailed diagnostics and correction, changes in the mean 
values and errors by the reduction of parameter values can be re- 
covered and considered for the interpretation of the code results. 
Additional uncorrectable deviations are only caused if the true 
error is smaller than the difference between adjacent parameter 
values. 

In summary, diagnostic plots with similar curves as pre- 
sented in this section can be used to evaluate whether mean val- 
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Fig. 11. Normalised diagnostic plot for the interpretation of 
mean values and standard deviations calculated by CIGALE 
for 10 equally-spaced parameter values between and 2, i.e. 
a step size of 0.22. The outer solid/dashed curves represent 
the enveloping curves for a Gaussian PDF having the mean in- 
side/outside the parameter range. The dash-dotted curve marks 
the enveloping curve for a box-shaped distribution. The triangu- 
lar area marked by "certain" indicates for a Gaussian PDF the 
region of negligible difference between apparent and real mean 
values (< 1% of the half parameter range) and errors (< 10%). 
For a box-shaped PDF this area would extend to higher cr and be 
delimited by the corresponding enveloping curve. For applying 
the diagram to real data, the right part of the parameter range has 
to be complemented by the mirrored diagnostic curves and the 
axes have to be multiplied by half the difference of the maximum 
and minimum parameter values selected. 



ues and errors resulting from CIGALE are reliable, uncertain, 
or unreliable (implying an unconstrained parameter). Fig. [TT] il- 
lustrates these different areas in a diagnostic plot similar to those 
used for the interpretation of the code results in Sect. |3.3.T] Only 
for data points in or close to the "certain" area the real PDF of a 
parameter can be well described by expectation value and stan- 
dard deviation provided by CIGALE, otherwise the shape of the 
PDF has to be studied in more detail or it has to be accepted that 
the parameter cannot be constrained for the photometric data and 
model set given. 



3. A test sample: SINGS 

In order to establish CIGALE as a tool for studying properties 
of nearby and distant star-forming galaxies, it has to be tested 
by using a well-known reference data set. In the context of the 
Spitzer Infrared Nearby Galaxy Survey (SINGS; Kennicutt et 
al. 2003) high-quality photometric data were obtained in the IR 
regime for a sample of representative nearby galaxies. Together 
with the available photometric data for the other wavelength 
ranges (Dale et al. 2007; Munoz-Mateos et al. 2009) this data 
set allows us to investigate the properties of the code in detail. 
We describe the sample in Sect. 13.11 and discuss the analysis of 
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Table 3. Selected model parameter values for the analysis of our 
SINGS sample 



10 20 
Distance [Mpc] 

Fig. 12. Optical (i.e. mainly R) absolute magnitudes and dis- 
tances in Mpc for the SINGS galaxies (see Kennicutt et al. l2003l l. 
The analysed objects are indicated by filled symbols. 



the data set by means of CIGALE in Sect. 13.21 The results and a 
comparison to the literature are presented in Sect. 13.31 

3.1. The sample 

The full SINGS sample consists of 75 representative nearby 
galaxies (Kennicutt et al. 2003), i.e., the different morpholog- 
ical classes from ellipticals to irregulars have similar weights. 
Moreover, the galaxies are distributed over roughly equal loga- 
rithmic bins in IR luminosity. Dale et al. (2007 2008 ) published 
UV-to-radio broad-bands SEDs for the SINGS galaxies. The 
Usted data comprises the FUV (~ 1500 A ) and NUV (~ 2300 A) 
filters of GALEX (Gil de Pa z et al. |2007] l, 2MASS data for J, H, 
and (JaiTett et al. l2003l l. and IRAC and MIPS data for the 
filters centr ed on 3.6, 4.5, 5.8, 8.0, 24, 70, and 160;um (Dale 
et al. I2005I I. Moreover, Dale et al. provide previously unpub- 
lished optical photometry for B, V, R, and /. However, since 
the optical data are erroneoufl we take the improved photome- 
try of Munoz-Mateos et al. ( I2009I I, which is based on asymp- 
totic magnitudes that were derived from surface photometry 
by extrapolating the curve of growth to encompass the whole 
galaxy light. Apart from all non-optical Dale et al. filters up to 
MIPS 160//m, their catalogue provides photometry of the Sloan 
Digital Sky Survey (SDSS; Stoughton et al. '2002') in u', g', r' 
i', and z' for 32 SINGS objects. Moreover, the Dale et al. (|2U07l 
120081 ) optical magnitudes of 23 galaxies without SDSS photom- 
etry were recalibrated by means of the catalogue of Prugniel 
& Heraudeau (119981 ), i.e., their aperture photometry was repli- 
cated on the SINGS optical images to derive zero-point correc- 
tions for each filter and each galaxy. For our sample selection we 
have required either the presence of SDSS photometry or at least 
three corrected Dale et al. filters. The latter criterion fulfil 16 of 
23 objects with recovered photometry only. Three of them do 
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not have 7?-band photometry. Finally, we only take galaxies for 
which all UV and IR filters are availabl^El Our selection criteria 
aim at ensuring a comprehensive photometric sample of similar 
high quality and wavelength coverage. Considering all criteria, 
our final sample comprises 39 SINGS galaxies (see Table |2]). 
Fig. [12] shows the basic properties distance and absolute mag- 
nitude (mainly R) as given by Kennicutt et al. (2003) for our 
subsample compared to the full sample. The diagram indicates 
that dwarf galaxies are underrepresented in our sample. There 
are no optical magnitudes for these galaxies in the catalogue of 
Munoz-Mateos et al. (120091 ) because of recalibration problems. 



3.2. Analysis 

In the following we discuss the model grid that was selected 
to analyse our test sample of SINGS galaxies (Sect. 13.2. It . 
Moreover, the quality of the fitting (Sect. 13.2.21 ) and the influ- 
ence of the filter set on the results (Sect. l3.2."3T ) is studied. 



3.2.1 . The model grid 

Studying the properties of our sample of 39 SINGS galaxies by 
SED fitting requires a careful selection of model parameters, 
since the parameter space has to be restricted for computing time 
reasons. Table [3] shows the model grid that we have eventually 
used. The individual values were selected to provide (roughly) 
equally-sized steps in the given units or in dex for dynamical 
ranges larger than one order of magnitude. Only eight param- 
eters have multiple input values, which is distinctly lower than 
the 15 to 17 filters available for the SINGS galaxies in our sam- 
ple. The total number of models is about 7 x 10^, which can be 
managed quite easily on a typical PC because of a runtime of 
CIGALE of a few hours only. 

For all stellar population models we use Maraston (120051) 
SEDs with Salpeter IMF and solar metallicity. The latter is jus- 
tified by a mean oxygen abundance (12 + log(0/H)) - 8.69 and 
a scatter of 0.16 for 23 sample galaxies listed in Calzetti et al. 



' The corrected photometry provided in Dale et al. ( 12008 1 is better 
than the original one given in Dale et al. ( 12007b . However, the results 
are still unsatisfying. 



Despite of the low quality of the far-IR data, two elliptical galax- 
ies were selected. NGC0584 was only marginally detected at 160yum. 
Its flux at this wavelength could be contaminated by diffuse back- 
ground/foreground emission. The fluxes at 70 and 160 fim of NGC 1404 
could be affected by a possible background source, which could not be 
masked and cleaned in a satisfying way due to the low image resolution 
at these wavelengths. 
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Table 2. Properties of the SINGS test sample. Although we list the relatively conservative run B results (see Sect. |3.2.2] i. the given 
standard deviations could underestimate the real uncertainties because of the limited model grid used, model inaccuracies, and 
possible unconsidered photometric errors. 
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±0.36 


NGC 3184 


SABcd 


8.6 


10.28 ±0.13 


0.11 


± 0.15 


-0.36 ± 0.06 


10.26 


±0.01 


9.76 


± 


0.04 


1.13 


±0.19 


NGC 3 190 


SAap 


17.4 


10.92 ± 0.03 


-0.85 


± 0.51 


0.82 ±0.12 


10.46 


± 0.01 


9.67 


± 


0.05 


2.06 


± 0.87 


NGC 3198 


SBc 


9.8 


10.10 ±0.12 


-0.03 


± 0.15 


-0.34 ± 0.04 


10.13 


±0.01 


9.60 


± 


0.04 


0.89 


± 0.17 


NGC 3351 


SBb 


9.3 


10.66 ± 0.08 


-0.04 


± 0.19 


0.17 ±0.18 


10.39 


±0.01 


9.83 


± 


0.05 


1.31 


±0.36 


NGC 3521 


SAB be 


9.0 


11.00 ±0.04 


0.43 


± 0.22 


-0.28 ± 0.03 


10.80 


± 0.01 


10.37 


± 


0.04 


2.74 


±0.37 


NGC 3621 


SAd 


6.2 


10.21 ±0.14 


0.12 


± 0.15 


-0.40 ± 0.09 


10.28 


± 0.02 


9.88 


± 


0.04 


1.31 


±0.26 


NGC 3627 


SABb 


8.9 


10.95 ± 0.05 


0.55 


± 0.22 


-0.29 ± 0.03 


10.79 


± 0.02 


10.40 


± 


0.05 


2.31 


±0.31 


NGC 4536 


SAB be 


25.0 


11.00 ±0.14 


0.96 


± 0.15 


-0.41 ± 0.07 


11.12 


± 0.02 


10.82 


± 


0.04 


2.00 


±0.33 


NGC 4559 


SABcd 


11.6 


10.27 ± 0.08 


0.12 


± 0.08 


-0.51 ±0.01 


10.45 


±0.01 


9.90 


± 


0.04 


0.80 


±0.18 


NGC 4569 


SABab 


20.0 


11.41 ±0.02 


0.49 


± 0.22 


0.33 ± 0.08 


11.05 


± 0.01 


10.29 


± 


0.07 


2.23 


± 0.37 


NGC 4579 


SABb 


20.0 


11.44 ±0.03 


0.23 


± 0.24 


0.44 ± 0.20 


11.04 


±0.01 


10.21 


± 


0.05 


2.02 


± 0.57 


NGC 4594 


SAa 


13.7 


11.75 ±0.02 


-0.47 


± 0.47 


0.92 ± 0.05 


11.27 


±0.01 


9.79 


± 


0.06 


1.60 


±0.89 


NGC 4625 


SABmp 


9.5 


9.35 ±0.13 


-0.72 


± 0.11 


-0.34 ± 0.05 


9.35 


± 0.02 


8.83 


± 


0.04 


0.86 


± 0.19 


NGC 4631 


SBd 


9.0 


10.47 ±0.14 


0.70 


±0.20 


-0.62 ± 0.08 


10.76 


±0.02 


10.50 


± 


0.04 


1.49 


± 0.21 


NGC 4725 


SABab 


17.1 


11.40 ±0.04 


0.35 


±0.18 


0.43 ±0.18 


11.01 


± 0.00 


10.14 


± 


0.04 


1.12 


± 0.33 


NGC 4736 


SAab 


5.3 


10.80 ± 0.07 


0.07 


± 0.20 


0.11 ± 0.23 


10.52 


±0.00 


9.89 


± 


0.04 


1.54 


± 0.36 


NGC 4826 


SAab 


5.6 


10.77 ± 0.02 


-0.32 


± 0.21 


0.31 ±0.10 


10.40 


±0.00 


9.60 


± 


0.04 


2.38 


±0.34 


NGC 5033 


SAc 


13.3 


10.77 ±0.11 


0.31 


± 0.20 


-0.27 ± 0.06 


10.67 


± 0.01 


10.30 


± 


0.04 


1.88 


±0.29 


NGC 5055 


SAbe 


8.2 


10.96 ± 0.06 


0.43 


± 0.22 


-0.28 ± 0.02 


10.78 


± 0.01 


10.33 


± 


0.04 


2.22 


±0.31 


NGC 5 194 


SABbe 


8.2 


10.82 ±0.12 


0.79 


±0.15 


-0.44 ± 0.07 


10.97 


±0.02 


10.62 


± 


0.05 


1.78 


±0.29 


NGC 5195 


SBOp 


8.2 


10.77 ± 0.02 


-0.12 


± 0.21 


0.33 ± 0.08 


10.41 


±0.01 


9.72 


± 


0.05 


3.36 


±0.38 


NGC 5474 


SAcd 


6.9 


9.30 ± 0.06 


-0.65 


± 0.14 


-0.45 ± 0.04 


9.50 


±0.02 


8.59 


± 


0.05 


0.27 


±0.06 


NGC 5713 


SABbep 


26.6 


10.77 ± 0.08 


0.73 


±0.19 


-0.53 ± 0.04 


10.93 


± 0.01 


10.71 


± 


0.03 


2.86 


±0.33 


NGC 5866 


SO 


12.5 


10.88 ± 0.03 


-0.92 


±0.54 


0.78 ±0.12 


10.43 


±0.00 


9.30 


± 


0.06 


2.51 


± 1.18 


NGC 7331 


SAb 


15.7 


11.39 ±0.08 


0.78 


± 0.20 


-0.04 ± 0.25 


11.17 


±0.02 


10.76 


± 


0.04 


2.66 


±0.50 


: Kennicutt et al. (|2003t 


























;nificant deviations of the metallicity from the solar value may affect the results for / 


D4000 (see Seet.|3.2.l|l. 











(I2007I I. Moreover, our sample lacks the low-metallicity dwarf 
galaxies present in the complete SINGS sample (see Fig. [T2] |. 
However, metallicity measurements are relatively uncertain due 
to high-quality requirements of the observations, insufficient 
coverage of a galaxy in combination with abundance gradients, 
and, in particular, calibration problems and general limitations 
of the metallicity tracer (see Moustakas & Kennicutt 2006l and 
references therein). Hence, we have checked the influence of the 
use of half and double solar metallicities on the results of the 
analysis of our SINGS sample. In general, the resulting prop- 
erties indicate deviations within the errors and the correspond- 
ing uncertainties increase by to 20% only if metallicity uncer- 
tainties are considered. The only exception is the effective age 
measured at 4000 A fD400o^ which indicates a decrease of about 
0.4 dex for a metallicity increase of a factor of 4. Consequently, 
the absolute values of fD40oo have to be taken with care due 



to the well-known age-metalUcity degeneracy (e.g., Kodama & 
Arimoto [T997t . 

In CIGALE SFHs are modelled by the combination of two 
models with exponentially decreasing SFR (see Sect. |2.1.TT i. For 
the "old" stellar population model we take a fixed age of lOGyr, 
which should be a good compromise between the age of the 
Universe and the most important phase of star formation for the 
galaxies investigated. The t values cover a relatively large range 
of values in order to account for the different kinds of SFHs pos- 
sible. Due to the low effect of old stellar populations on galaxy 
SEDs, it is sufficient to select a few values only. For the "young" 
stellar population model we also take a few ages and r only, 
since the variety of SEDs dominated by young stars is not very 
large. On the other hand, the mass fractions of both t models are 
critical. Therefore, we analyse nine burst fractions /burst covering 
the full range of theoretically possible values. 
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^GLX 

Fig. 13. The IR-to-far-UV luminosity ratio Lir/Lfuv versus the 
UV continuum slope yScLX for the SINGS sample. The solid 
curve indicates the theoretical position of starbursts in the dia- 
gram (Kong et al. "2004'). The area in the yScLX-iiR/ipuv plane 
that is covered by CIGALE models is marked by small dots. 
Subfigure (a) shows the full model set for the parameters given 
in Table [3] In (b) no slope correction of the Calzetti law is con- 
sidered, i.e. 6-0. Finally, in (c) the stellar radiation is reddened 
regardless of the stellar population, i.e. /att - 1. 



The obscuration of the stellar populations by dust is consid- 
ered by the application of attenuation laws with different reason- 



able slopes (i.e. slopes not too far from the Calzetti et al. ( 120001 ) 
law) but without a UV bump. Since the UV wavelength range of 
the SINGS galaxies is covered by the two GALEX filters only, 
details of the attenuation law are difficult to study and insignifi- 
cant for most galaxy properties as test runs of CIGALE indicate 
for a wide range of attenuation-related parameters. Therefore, 
we can restrict the number of models by only taking the default 
value zero of the Ry-ielaled parameter cr and by only consider- 
ing slope corrections 6 between -0.3 and 0.3, which result in real 
Rv between 3.0 and 5.9 (see Sect. 12. 1.2] and Fig.[3]l. Moreover, 
we fix the UV bump strength and select a value of zero. Apart 
from the negligible effect of the UV bump on the fit results for 
a wide range of strengths (cf. Noll et al. ,2009) , this choice is 
justified by the non-detection of a UV bump in local starburst 
galaxies (Calzetti et al. 1994), the lack of a significant 2175 A 
feature in the Kinney et al. (1996) characteristic UV spectra of 
nearby galaxies covering a large range of star formation activity, 
the only low-to-moderate UV bump strengths in different galaxy 
populations (also comprising "normal" star-forming galaxies) 
studied at intermediate/high redshift (Noll et al. 2009), and the 
predicted weakening of the UV bump in local spiral galaxies by 
radiative transfer effects (Silva et al. |1998t Pierini et al. 120041 ). 

Since the V-band attenuation of the young r model Ay ^sp 
has an important effect on the (UV) SED slope and the UV-to- 
IR flux ratio, we choose a fine grid of values and a sufficiently 
high upper limit. The attenuation factor /an indicates the fraction 
of Av.ysp valid for the old t model. The possibility to vary this 
parameter is crucial for fitting non-starburst galaxies as typically 
found in our SINGS sample. In Fig. [13] we show the luminosity 
ratio Lir/Lfuv and the UV continuum slope /3clx of our sample 
galaxies. Lir was derived from the MIPS 24, 70, and 160 jum fil- 
ters as described in Dale & Helou ( |2002| l. And based on GALEX 
FUV and GALEX FUV - NUV, respectively, Lpuv and /3glx 
were calculated following the recipes of Kong et al. (2004). The 
objects are located close to or below the so-called "Meurer" re- 
lation (Meurer et al. |1999t Kong et al. 2()04). This curve is valid 
for starburst galaxies where the stellar content is completely ob- 
scured by a dust screen or shell. As the results of CIGALE in 
subfigure (c) imply, the galaxy distribution found cannot be ex- 
plained by a deviation of the effective attenuation curve from 
the Calzetti law. Instead, subfigure (b) illustrates that an age- 
dependent amount of attenuation (cf. Silva et al. |1998t Pierini 
et al. 2004; Panuzzo et al. 2007; Nofl et al. 2007) as simulated 
by /att is able to reproduce the distribution of SINGS galaxies. 
Then, the data points that deviate most from the starburst curve 
are characterised by high obscuration of a young stellar compo- 
nent of very low mass (low /burst) and nearly unattenuated old 
stellar populations. Consequently, our code is suitable to study 
"normal" star-forming galaxies which are usually located below 
the Meurer curve and are typical of the nearby Universe (see 
B uat et a l. 120051 Cortese et al. l2006l Dale et al. l2007t Salim et 
al. lIOOTl l. 

Dale & Helou ( 120021 ) offer IR dust emission templates in 
the a range from 0.0625 to 4.0. Our selection covers a values 
between 1.0 and 3.0. The narrower range is reasonable, since 
cc > 1.0 correspond to mean IR luminosities Ljr ^ 4 x lO'-' 
according to the calibrations of Chapman et al. (2003) and 
Marcillac et al. ( I2006I I. This is a safe limit even if the scatter in 
the calibrations is considered, since the SINGS sample is char- 
acterised by relatively low IR luminosities (see Kennicutt et al. 
120031 ). For very high a values the IR SEDs become very similar 
for the wavelength range studied. The calibrations of Chapman 
et al. and Marcillac et al., which are based on the ratio of the 
rest-frame fluxes at 60 and lOOjum, are degenerated in this case. 
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Table 4. Mean photometric errors for our SINGS test sample 



Filters 


Rel. errors 


r^AT fIJV \J1J\/ 


IJ/C 


Dale et al. B, V, R, I 


16% 


SDSS u', g', r', i', z' 


3% 


2MASS J, H, K, 


1% 


IRAC3.6, 4.5,5.8, S.Ofim 


11% 


MIPS 24fim 


5% 


MIPS 70 ^m 


7% 


MIPS 160/im 


13% 



0.6 



^ 0.4 
to 



0.2 



o. 



-0.4 



~1 T 1 T 

Run A 



"1 r- 



"I r 



"T 1 1 r 



_L 



_L 



-10 12 
log X [yLtm] 

Fig. 14. Ratio of the best-fit model fluxes and the measured ob- 
ject fluxes in mag for run A and different filters and objects. The 
filters are indicated by their mean wavelength. Sample-related 
mean and median deviations for each filter are marked by big 
circles and lozenges, respectively. The open symbols indicate 
the recalibrated optical filters of Dale et al. (I2007i) . 



Therefore, it does not make sense to analyse a values beyond 
3.0. 

Finally, we do not use the option of an AGN contamination 
or hot-dust contribution for the main run. Although there are sev- 
eral Seyfert galaxies in the sample (Kennicutt et al. '2003'), we 
do not expect that the total magnitudes of the selected galaxies 
are significantly affected by such a contribution. Nevertheless, 
we have checked with a different set of models allowing for a 
large range of hot-dust IR contributions whether CIGALE repro- 
duces the expected result. Indeed, the vast majority of IR SEDs 
is consistent with no contribution at all. Only for the early-type 
galaxy NGC 1404 we find a significant (> 10%) hot-dust frac- 
tion. However, this result could also be explained by possible 
difficulties to fit the mixture of stellar radiation and very weak 
dust emission in the mid-IR of this special galaxy. Moreover, 
the fluxes at 70 and 160//m have to be taken with care due to 
a possible contamination by a background source (see note in 
Tabled. 




-4 
2 



NGC4625 




J I I I I I I L 



4 6 
log X [A] 

Fig. 15. Illustration of three characteristic best-fit SEDs (run A) 
which differ in their UV-to-optical ffux ratios. From top to 
bottom, the SEDs of NGC 4594 (Sa), NGC 5033 (Sc), and 
NGC 4625 (Sm) are shown. The measured photometric ffuxes 
are marked by circles. 



3.2.2. The fit quality 

For the analysis of the properties of our SINGS test sample we 
use the total magnitudes of 21 filters and its errors reported in 
Munoz-Mateos et al. ([2009). The typical uncertainties for the 
different filters are listed in Table |4] They are dominated by the 
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instrument-related zero-point errors or uncertainties of photo- 
metric corrections and vary between 1% for the 2MASS near- 
IR filters and 16% for the Dale et al. (2007) optical magnitudes 
corrected by Munoz-Mateos et al. ( I2009I I. The latter are only rel- 
evant for 15 of 39 sample galaxies for which SDSS photometry 
of much higher precision is not available. Since the relative pho- 
tometric errors represent the weights for the different filters in 
the fitting process, the weights of the optical magnitudes for the 
two subsamples with and without SDSS photometry are quite 
different. Therefore, we will discuss the results of both subsam- 
ples separately in Sect. 13.31 

The large differences in the reported uncertainties for the fil- 
ter fluxes obtained by different projects and instruments raises 
the question whether these values are really comparable. The 
homogeneity of the error evaluation is crucial for reliable fit- 
ting results. In particular, very low errors such as reported for 
the 2MASS filters are critical, since this causes hard constraints 
for the model set. Moreover, especially small errors are ques- 
tionable, since possible systematic offsets between different fil- 
ters can cause erroneous fits. Since the models are not a perfect 
reproduction of the nature, systematic errors are also expected 
for them if the fluxes of different wavelength ranges are com- 
pared. For most wavelengths we assume uncertainties in the or- 
der of 5 to 10%. Higher values are especially expected for the 
rest-frame wavelength range between 2.5 and 5 fim for which 
the continuum is relatively uncertain due to the mixing of stel- 
lar and dust emission and fundamental uncertainties in the shape 
of both components. An example is the insufficient knowledge 
of the contribution of TP-AGB stars at these wavelengths (see 
Maraston 2005). Therefore, we study two different photometric 
catalogues. In the first case ('run A') we take the published er- 
rors (see Table m and in the second case ('run B') we increase 
them by adding an additional moderate 5% error in quadrature 
to all filter errors. This procedure allows us to investigate how 
the results change if additional systematic uncertainties are con- 
sidered. 

As a first result of CIGALE, Fig.[T4]shows the deviations of 
the best-fit models from the object photometry in magnitudes for 
each filter Since the best-fit models of both runs indicate similar 
photometric fluxes, we only plot the data of run A. In general, the 
sample-averaged deviations for the different filters are relatively 
small. On average, the difference is 0.07 mag only. The worst 
filters are MIPS 160yum, IRAC 3.6 yum, and IRAC 4.5 yum with 
deviations amounting to 0.28, 0.19, and 0.17 mag. The modest 
fit quality for these three filters is caused by relatively high un- 
certainties in the object photometry and the shape of the models 
in the corresponding wavelength ranges (see above). For MIPS 
160/im the significant deviations can partly be explained by the 
lacking flexibility of the one-parameter models of Dale & Helou 
( I2002I I in reproducing the complex IR SEDs of real galaxies. 
On the other hand, the differences between the Dale & Helou 
templates of adjacent a values (see TableO are distinctly larger 
for MIPS 160/im than for any other IR filter studied. The corre- 
sponding uncertainties are similar to those indicated in Fig.[T4l 
This error source becomes negligible if expectation values de- 
rived from the parameter PDFs are analysed instead of best-fit 
models. 

The good fit quality that we can reach with the parameter 
set chosen over a wide wavelength range is demonstrated in 
Fig. [15] which shows the best-fit SEDs of three SINGS galaxies 
in comparison to the measured photometry. Although the dust 
reddening and star formation activity of NGC4594, NGC 5033, 
and NGC 4625 differ considerably, the photometric SEDs of all 
galaxies shown are reproduced quite well. 




_Q 2 I I I I I I I I I I I I I I 

10.5 10.6 10.7 

log Mstar [%] 

Fig. 16. Effect of the reduction of IR filters on the basic galaxy 
properties total stellar mass and SFR. The diagram provides the 
sample means and mean errors of these quantities for run A 
(solid lines) and run B (dashed lines) and different filter com- 
binations. The labels near the symbols indicate the mean wave- 
length of the last filter in the filter set in fim. 160 means that 
all filters from GALEX FUV to MIPS 160 jum are used for the 
analysis. 2.2 refers to the band as the reddest filter 



Finally, we have checked whether our approach to derive 
PDFs has any influence on the results. In general, we find 
that the differences between the 'max' and 'sum' methods (see 
Sect. 12.2b are negligible, i.e., the differences in the expecta- 
tion values are distinctly smaller than the errors. For example, 
the 'sum' method (for a lower model confidence limit of 10"^) 
changes the sample-averaged total stellar mass and SFR by only 
-0.01 dex (cr = 0.03) and +0.04 dex (cr = 0.09) for run A, re- 
spectively. In comparison, the corresponding errors are 0.05 dex 
and 0. 16 dex. They increase by 12% and 2%, respectively, by the 
use of the 'sum' method instead of our preferred 'max' method. 
Although the sample-averaged differences are small, significant 
effects cannot be excluded for individual galaxies, however The 
most extreme deviation for the SFR is an increase of 0.34 dex for 
NGC 5866. It is caused by diminishing, additional peaks in the 
probability distribution at low SFRs due to a low number of mod- 
els populating these peaks. Although such cases are rare as the 
comparison of the galaxy properties derived by both approaches 
show^3> it is the reason why we prefer the 'max' method, which 
reduces the influence of the model density on the weighing of 
parameter values. 

3.2.3. Influence of filter set on results 

Since CIGALE allows us to fit photometric data of galaxies rang- 
ing from far-UV to far-IR, star-formation-related parameters can 
be derived in a reliable way due to a full coverage of the dust- 
affected wavelength range. However, for distant galaxies espe- 
cially far-IR data are often not available or not sensitive enough. 



" The differences in the SFRs of the 'max' and 'sum' methods are 
relatively large in comparison to the results for other model parameters. 
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Table 5. Change of the SINGS sample mean for different prop- 
erties by running CIGALE without IR data 



Par. 


Run A 


Run B 


logMst^r 


n 1 n 


-U./O 


log SFR 


+0.83 


+0.67 


log Lboi 


+0.37 


+0.31 


log Ldust 


+0.75 


+0.15 


log /burst 


+ 1.05 


+0.72 


log ;d4000 


-0.78 


-0.71 


^FUV 


+2.16 


+ 1.59 


Av 


+0.65 


+0.57 


a 


-0.40 


-0.44 



Therefore, we have investigated how results would change if the 
IR filter set of our sample galaxies was incomplete. 

Fig. [16] shows the sample-averaged total stellar mass and 
SFR for a reduction of the IR filter set. While the full filter set 
contains 160jum as the reddest filter, the most reduced filter set 
ends at 2.2//m, i.e. the K^^ band. There is no significant change 
of the results as long as the filter set still includes at least one fil- 
ter (i.e. 24yum) which traces the dust emission beyond the PAH 
regime. Amplitude and slope of the Dale & Helou (2002) mod- 
els appear to be well determined in this case. In contrast, the 
sample means change considerably if the IR is only covered up 
to 8 fim, i.e. the region of the PAH emission features, where the 
Dale & Helou models show little variety only. For run A, Mstm 
decreases by 0.19 dex and the SFR increases by 0.83 dex if K^^ 
is the reddest filter These numbers show that the disregard of 
mid-IR and far-IR data results in unsatisfying fits for our SINGS 
test sample. Table|5]shows that many parameters are concerned. 
Apart from the SFR, important changes are also present for the 
dust luminosity (+0.75 dex), the burst fraction (+1.05 dex), and 
fD4000 (-0.78 dex). The effective attenuation factors in the far- 
UV Apuv and visual Ay considerably increase by 2. 16 mag and 
0.65 mag, respectively (cf. Burgarella et al. |2005l l, if IR data are 
not considered. Compared to other properties the effect on the 
mass is relatively small, since this quantity is only indirectly af- 
fected by a higher but still minor importance of young stellar 
populations and the corresponding change of the mass-to-light 
ratio. 

Fig- m and Table |5] also illustrate the results for higher un- 
certainties in the filter fluxes (run B), which indicate lower de- 
viations for most properties. In particular, the changes for Ljust 
(+0.15 dex), /buret (+0.72dex), andApuv (-i-1.59mag) are signif- 
icantly reduced. This suggests that the low photometric errors in 
some filters of run A could force the code into unfavourable fits 
if no IR data is available. Consistently, the galaxies with SDSS 
data and consequently lower photometric errors in the optical in- 
dicate higher differences in the mean SFR of code runs with and 
without IR data than objects with corrected optical Dale et al. 
photometry only. 

We have also studied the code results for a filter set with- 
out the GALEX FUV and NUV filters. The differences for most 
parameters are insignificant and the SFR increases by 0.17 dex 
(run A) or 0.11 dex (run B) only. Even Apuv that directly de- 
pends on the UV flux changes by +0.02 mag for run A only. An 
exception is the increase of this quantity by 0.46 mag for run B. 
However, this amount is still much lower than the deviations 
found for the neglection of IR data. Consequently, the UV data 
of the SINGS galaxies have only little influence on the fit results. 
This is probably caused by the relatively low weight of these fil- 



ters in the fitting process (see Table |4]i. The relatively high un- 
certainties in the photometry of these two filters do not appear to 
significantly constrain the star-formation-sensitive UV SEDs of 
the sample galaxies. Consequently, a lack of IR data can cause 
striking systematic deviations in the galaxy properties. However, 
the amount of these deviations probably depends on the SFHs of 
the galaxies investigated and the filter set available. For example, 
missing IR data could be better compensated at higher redshifts 
because of a better filter coverage in the rest-frame UV due to 
the shift of the observed frame of optical telescopes. 

3.3. Results 

In the following we present the physical properties of our SINGS 
sample obtained by means of CIGALE and discuss their reliabil- 
ity (Sect. 13.3. 11 1. We also show relations between the different 
parameters and compare our results to those from other studies 
(Sect. Ilia . 

3.3.1. Expectation values and standard deviations 

In Fig. [17] (run A) and Fig. [18] (run B) we illustrate the expec- 
tation values and standard deviations for the different model pa- 
rameters of CIGALE for our test sample of 39 SINGS galaxies. 
We only show results for parameters with PDFs consisting of 
at least two values. Data points based on photometry with and 
without SDSS photometry are marked by different symbols (cir- 
cles and crosses). For evaluating the reliability of the results, the 
diagrams with mass-independent properties also indicate the di- 
agnostic curves introduced in Sect. 12.31 In fact, the enveloping 
curves for Gaussian and box-shaped PDFs as well as the "safe" 
areas for Gaussian distributions are shown. 

At first, we discuss the mass-dependent properties Mstar, 
SFR, Lboi, and Ldust (see Table |2|. Typical total stellar masses 
of about lO'^Mo and SFRs between 0.1 and lOMoyr^' re- 
flect that the sample mainly consists of nearby spiral galaxies 
(see Sect. 13.1b . Moreover, the dust luminosities are lower than 
10" Lq, which indicates that luminous IR galaxies (LIRGs) are 
not present in our sample. The mass-dependent properties are 
quite well constrained (see Fig.|6]for an example PDF). The typ- 
ical run A errors range from 0.011 dex for Lboi to 0.16 dex for 
the SFR, which is distinctly smaller than the range of parame- 
ter values covered by the sample. This statement is also true for 
the somewhat higher average errors between 0.013 and 0.22 dex 
for the less constrained run B. The differences between the pa- 
rameter values for run A and run B are usually within the errors. 
For the most uncertain mass-dependent property, the SFR, the 
sample-averaged value increases by 0.05 dex only if the SFRs of 
run B instead of run A are considered (see Fig. [T6] l. The mean 
difference between the SFRs of both runs for individual galaxies 
is 0.10 dex. 

In contrast to basic model parameters which describe the 
SFH or the attenuation of the stellar radiation (see Table [T), the 
mass-dependent global galaxy properties are relatively indepen- 
dent of the details of the models used. Hence, the results for our 
SINGS test sample can be compared quite comfortably to other 
studies in this cas43 The dust luminosity Ldust is relatively easy 
to measure, since it is directly related to the filter fluxes in the 
IR (provided that the filters cover the entire relevant wavelength 



Unfortunately, the completely different parameter set and a signifi- 
cantly different selection of test galaxies from SINGS impedes a qual- 
itative and quantitative comparison of our code results to those of da 
Cunha et al. I i2008i ). 
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Fig. 17. Run A code results for our SINGS test samples with and without SDSS photometry (circles and crosses, respectively). 
Expectation values and standard deviations are shown for mass-dependent properties (Mstar, SFR, Lboh and Ljust), basic model 
parameters with at least two selected input values (tqsp, Tysp, fysp, /burst, S, Ay^ysp, /att, and a), complex ages (tu and fD4()Oo), and 
effective dust attenuation factors (Apuv and Ay). For the basic parameters the selected parameter values are indicated by vertical 
dotted lines. The different diagnostic curves are described in Fig.[TT] All curves depend on the number of parameter values given. 
For the complex ages and effective attenuation factors we assume 50 which is the number of bins of the parameter probability 
distributions. The curves for the mass-dependent parameters are not shown, since the errors are much lower than the widths of the 
object-dependent parameter ranges. 



range). Nevertheless, the comparison of Ldust obtained by differ- 
ent studies is a good consistency check. For 13 sample galaxies 
Ldust determinations are also available by Draine et al. ( I2007I I 
who derived them from the SINGS IR SEDs using the Draine 
& Li ( I2007I I dust models. The mean difference between our es- 
timates and those of Draine et al. is negligible, since it results 
in 0.00 ± 0.03 dex for run A and B. The scatter amounts to 
0.1 dex. A comparison of instantaneous SFRs is more crucial, 
in particular, if the approaches are completely different. For 29 
sample galaxies SFRs of Kennicutt et al. (120031) are available. 
They are based on Ha emission, which is a thoroughly studied 



and relatively reliable SFR indicator (see, e.g., Kennicutt |1998t 
Brinchmann et al. |2004l Kewley et al. |200?l Calzetti et al. |2007] l. 
Despite of the completely different methods, our SFRs of run B 
exhibit a moderate relative scatter of 0.3 dex and deviate on av- 
erage by 0.06 + 0.05 dex only from those of Kennicutt et al. (see 
Fig. [T9] ). For run A the SFRs agree even better, since the mean 
difference amounts to 0.00 + 0.06 dex. In any case, this compari- 
son shows that our results for our SINGS test sample are reliable 
at least for galaxy properties which are relatively independent of 
the details of the modelling. 
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Fig. 18. Alternative run B code results for our SINGS test samples. In comparison to Fig. [17] the given errors of the input photometry 
were increased by 5% added in quadrature in order to account for additional systematic photometric errors and model-inherent 
uncertainties. 



The basic stellar-population-related parameters are the ages 
and T of two Maraston ( 120051 ) models, which are linked by the 
burst fraction /burst (see Sect. 12.1. H and TableO. We show them 
in Figs.[T7]and[T8]excepting the constant f osp- For Tosp and run A 
the resulting values cover almost the entire investigated range 
between 0.25 and 10 Gyr. However, the errors are relatively high. 
Most of the data points are close to or on the enveloping curves. 
In particular, the distribution appears to follow the enveloping 
curve for box-shaped PDFs. If those distributions were typical 
(which is difficult to prove because of the small number of pa- 
rameter values), the Tosp of most galaxies could be classified as 
"rather low" or "rather high" only. However, the situation ap- 
pears to be even worse if the results of run B are considered. The 
lower accuracy of the input photometry causes that for the ma- 
jority of the objects Tosp is not constrained at all, as the cluster- 
ing of data points in the peak of the enveloping curves indicates. 



Consequently, it is safer to assume that the given photometry is 
not accurate enough to derive t values for the old stellar pop- 
ulation model. This also appears to be the case for the age and 
T of the young stellar population model. There is only a trend 
towards relatively high mean f ysp for run A. However, this trend 
vanishes completely if the results of run B are considered. In 
contrast, the burst fractions of most sample galaxies are reliable 
at least in the context of the SFHs allowed by the model grid. 
Even for run B the data points cluster inside the area of trust- 
worthy results. If run A provided the correct results, the errors in 
/burst of part of the galaxies would be even below the grid spac- 
ing, which is suggested by a position below the "triangle" in the 
diagnostic diagram. In any case, small burst fractions below 1 % 
dominate the distribution. No galaxy shows /burst > 10%. This 
result implies that the current star formation activity tends to be 
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log SFR (Kennicutt+03) [Mo/yr] 

Fig. 19. Comparison of the SFRs derived from run B of CIGALE 
for the SINGS sample and those of Kennicutt et al. (2003 ) based 
on Ho- emission. Upper limits are indicated by arrows. The error 
bars of the code results indicate the standard deviation of the 
SFR probability distribution for each galaxy. 



lower than the past average, since the age ratio of the old and 
young stellar population model is about 80. 

The relatively high degree of degeneracy in the SEDs for dif- 
ferent SFHs and the enormous amount of possible evolutionary 
scenarios prevents the reliable derivation of most basic model 
parameters related to the stellar population. Hence, we also dis- 
cuss the effective ages tu and f 04000 (see Sect. |2.1.TT l. The mean 
mass-weighted age Im of 8 Gyr is close to the age of the old stel- 
lar population model of 10 Gyr. Consequently, the galaxies un- 
der study appear to be dominated by old stars that were formed 
in the first half of the lifetime of the Universe. However, the ex- 
act fM of the individual galaxies is relatively uncertain. While 
one third of the objects indicates a trustworthy age in the case 
of run A, there are only a few galaxies with reliable ages in the 
case of run B. The decrease in errors around 6 Gyr is related to 
the maximum burst fraction of 10% found in the sample. The 
age fD4ooo derived from the strength of the 4000 A break indi- 
cates a wider variety than f m and is quite well constrained even 
for run B (see also Table |2]i and even under consideration of 
the metallicity dependency of the results (see Sect. l3.2.TI ). If the 
mean fD40oo of about 1 Gyr is taken, the sample can be divided 
into a spectroscopically young subsample of 22 galaxies and a 
mean age of 400 Myr and a complementary old one of 17 galax- 
ies with 3.6 Gyr on average. The relatively wide range in f d40oo 
reflects the variety of galaxy types in SINGS ranging from ellip- 
ticals to irregulars (see Kennicutt et al. 120031 ). While 90% of the 
young subsample indicates morphological types of Sb and later, 
only 20% of the old subsample exhibits these types. 

Since we do not scale the attenuation law by a Ry-iike con- 
stant factor and do not consider the possible presence of 2175 A 
absorption features in the SEDs (see Sect. 13.2. It , the shape of 
the attenuation curve is only determined by the deviation 6 of the 
slope from a Calzetti law. Figs. [TtI and [TS] show that 6 of many 
galaxies is not constrained at all and the rest is highly uncertain. 



Interestingly, significant deviations of 6 from zero are almost ex- 
clusively found for run A and galaxies with SDSS photometry. 
The latter has not played a role for the parameters discussed so 
far This result implies that differences in the shape of the atten- 
uation curve can only be detected if the accuracy of the photom- 
etry is in the order of a few percent. For the SINGS sample the 
accuracy is not good enough for any reliable statements. The sit- 
uation would probably be better if more and better UV data were 
available for the SINGS galaxies. However, this is not the case. 
In comparison to 6 the scaling parameter Ay_ysp for the young 
stellar population, which ranges from 0.15 to 1.66 for run A, 
can be constrained quite well. Most galaxies are in the reliable 
region of the diagnostic diagram. However, the individual val- 
ues significantly depend on the run, in particular, if SDSS pho- 
tometry is available. In this case, the mean Av,ysp changes from 
0.67 for run A to 0.87 for run B. For the other galaxies the aver- 
age increase is only half as strong (from 0.89 to 1.00). This be- 
haviour can be explained by the stronger run A SED constraints 
in the optical for galaxies with SDSS photometry and the cou- 
pling of 5 and Ay ysp that is suggested by the stability of the dust 
luminosity. The steeper slopes of the attenuation law for run A 
({6} = -0.04 versus -0.01 for run B) have to be compensated by 
lower Av ysp in order to reach a similar obscuration of the stel- 
lar populations by dust. Another attenuation-related parameter 
is /att which gives the fraction of Av,ysp for the old stellar pop- 
ulation. For individual galaxies the interpretation of the results 
is uncertain especially for run B. On the other hand, the mean 
/att for run A and run B are 0.45 + 0.03 and 0.41 + 0.02, respec- 
tively. This suggests a trend towards /an < 0.5 in SINGS. The re- 
sult is consistent with the fact that many galaxies in our SINGS 
sample are highly diversified galaxies for which the amount of 
dust attenuation depends on the age of the stellar population (cf. 
Fig. [nil. 

The results for the attenuation at 1500 A are reliable for most 
galaxies. The parameter space allowed by the model grid is only 
partly covered by sample objects. Relatively low Apuv dominate. 
The mean values are 1.5 (run A) and 1.6 (run B). They corre- 
spond to 78% and 69%, respectively, of the values expected for 
the sample Ay ysp and the Calzetti law (S - 0), which is charac- 
terised by Apuv/Ay - 2.55. Hence, the old r model (f/r ~ 10) 
significantly contributes even in the far-UV, which is obviously 
due to the very low burst fractions and low /att of most sample 
galaxies. The effective attenuation factors in the visual Ay are 
less reliable than Apuv and Ay ysp. In particular, the run B re- 
sults are relatively uncertain. This outcome can be explained by 
the strong dependence of Ay on the relatively uncertain atten- 
uation factor /att. The mean Ay of both runs is about 0.3 mag, 
which corresponds to mean Ay/Ay ysp of 0.38 (run A) and 0.29 
(run B). These values are comparable to the /att of the sample 
and suggest that the optical continuum of most sample galaxies 
is almost completely given by the old r model. 

Finally, we discuss the shape of the dust emission in the IR. 
Since we do not consider a hot dust continuum (see Sect. 13.2. It , 
it is only determined by the a parameter of the Dale & Helou 
(I2OO2I 1 templates (see Sect. |2.1.3l l. Figs.[T2]and[T8]indicate that 
low a are very well constrained, while high a are relatively un- 
certain. This is obviously caused by the degeneracy of the dust 
emission models for wavelengths shortwards of the emission 
peak for high a. For c > 2.5 the flux ratio /v(60 fj.m)/ fy(lQO jjm) 
becomes almost constant. Relatively high a dominate the dis- 
tribution, i.e., the dust temperature tends to be cool. The aver- 
age sample value of 2.4 can be compared to the mean Ljust of 
6 X IO^Lq by using the calibrations of Chapman et al. (I2003I I 
and Marcillac et al. (i2006j for which we obtain a x 2.1 and 
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Fig. 20. Relations between the total stellar mass Mstar, the specific SFR 0, the effective age at 4000 A fD4()00, and the attenuation in 
the far-UV 

^Fuv for the SINGS test sample. The given values and errors are related to run B of CIGALE. The different symbols 
indicate different morphologies (see legend in upper left subfigure). The solid and dashed curves in the plots not related to Mstar 
show the mean values and 1 cr scatter of the model distribution. 



2.3. The moderate differences in a could be explained by differ- 
ent properties of the samples used for the calibrations and the 
SINGS sample. Differences could also be caused by the degen- 
eracy in /v(60 yum)//v(100 /urn) for high a values and the fact that 
this flux ratio was taken for the Chapman et al. and Marcillac et 
al. calibrations instead of a. 



3.3.2. Correlations 



The different galaxy properties derived by CIGALE are not com- 
pletely independent of each other. For example, the burst fraction 
/burst and the two effective ages ?m and f d4000 are well correlated, 
since the stellar population properties of the models used mainly 
depend on the burst fraction (see Sect. 13.3.11 1. Another example 
is the relation between the slope modification of the Calzetti law 
6 and the attenuation in the visual of the young r model Ai/ysp 
based on a relatively stable dust luminosity Ldust discussed in 
Sect. 13.3.11 Finally, Mstar, SFR, Lboi, and Ljust tend to be cor- 
related because of their dependence on the galaxy mass (see 
Table [1]). Therefore, we restrict our discussion of possible cor- 
relations between different galaxy properties on parameters for 
which the models do not show a direct relation. Moreover, we 
only consider parameters which were well determined in run B 
of CIGALE for our SINGS sample (see Table |2]i and, therefore, 
do not suffer from crucial degeneracies. In fact, we discuss the 
CIGALE output parameters M^tar, fD4000, and Apuv- Moreover, 



we show the specific SFR 0, i.e. the instantaneous SFR divided 
by the total mass produced by star formation in the pasQ 

Figure |20] presents the mutual relations of Mstar, 0, fD4000, 
and Apuv for our sample of 39 SINGS galaxies. The tightest 
correlation is found for the quantities and fD4000, which de- 
pend on stellar population parameters only. It shows that effec- 
tively younger stellar populations are linked to higher star for- 
mation activity. The correlation obtained is in good qualitative 
and quantitative agreement with the results of Brinchmann et al. 
(120041) for the relation between Ha-based (p and D4000 of star- 
forming SDSS galaxies. However, our distribution is narrower, 
since a tight correlation is already inherent in the code and the 
model grid used as the distribution of models in the 0-fD4OOO 
plane indicates. For the Mstar-related diagrams in Fig.|20lcorre- 
lations based on the selected model grid or the internal structure 
of the models are not possible, since the mass is the free scaling 
parameter in the fitting process. Therefore, all mass-related cor- 
relations found must have an astrophysical interpretation only. 
Although the correlation is weaker than those discussed before, 
the specific SFR clearly tends to decrease with increasing mass, 
which is in qualitative agreement with previous studies (see, e.g., 
Brinchmann et al. 2004; SaUm et al. 2007; Buat et al. 20(21 
2008). Quantitatively, there are some differences due to different 
sample properties. For example, Salim et al. (120071 ) find an aver- 
age (p for their sample of star-forming SDSS galaxies at z «; 0.1 



'^^ For Maraston ( 12005 1 models the time-integrated SFR corresponds 
to the galaxy mass, which comprises the total stellar mass and the mass 
of gas released from stars by winds/explosions (see Sect. |2. l.TT ). 
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Fig. 21. Comparison of the attenuation in the far-UV Apuv 
(run B) and the luminosity ratio Lir/Lfuv as defined by Kong 
et al. (I2004I I for our SINGS test sample. For the meaning of the 
symbols see Fig. |20] Errors can only be shown for the code re- 
sults. The solid curve exhibits the relation between both parame- 
ters for a solar metallicity Maraston ( i2005j model with constant 
SFR and an age of lOOMyr that is reddened by a Calzetti law. 
The dashed curve shows the corresponding relation for nearby 
star-forming galaxies of Burgarella et al. (i2005j . 



that is about0.4dex higher at 10 Mq than the value of -1.3 dex 
typical of our sample. The transition from active star forma- 
tion to a relatively passive evolution is at about 10" Mq for our 
SINGS galaxies. This break in the galaxy properties is even more 
pronounced for the effective age f d4()00- The bimodality found is 
consistent with the results of Kauffmann et al. (I2003a[ I2003bl l 
for the SDSS and is reminiscent of the two basic galaxy pop- 
ulations "red sequence" (> 7 objects) and "blue cloud" (> 19 
objects) discussed in, e.g., Strateva et al. (1200 11 1. Baldry et al. 
( |200 4i and Driver et al. ('2006'). The galaxies in-between have 
the highest age uncertainties and could represent at least in part 
the so-called "green valley". 

In Fig.|20]we also identify the morphological types of our 
sample galaxies as given in Table|2](see Kennicutt et al. I2003I I. 
For the relations discussed so far, the locus of a galaxy is clearly 
related to the morphological properties. The earlier the type of 
a galaxy is, the lower (p and the higher fD40oo and Mjtar are. 
These correlations show that our code results are consistent with 
the well-known dependence of the star formation activity and 
mass of nearby galaxies on the morphological type (see Roberts 
& Haynes 1994 and Kennicutt 1998 for reviews). In this con- 
text, the extreme properties of the early-type SBa spiral galaxy 
NGC 2798 can be understood by the presence of a nuclear star- 
burst that produces a high amount of deeply dust-enshrouded 
young stars. 

No convincing correlation is found for the attenuation at 
1500 A Apuv and (p or fD4()00, i-C- the obscuration of the young 
stellar population is relatively independent of the star forma- 
tion activity and history. Concerning the morphology, the lat- 
est types (and possibly the ellipticals) in the sample tend to 
have lower Apuv than most spiral galaxies. For instance, Apuv 



equals 2.00 + 0.29 for SO-Sab, 1.64 + 0.17 for Sb-Scd, and 
0.92 + 0.22 for Sd-I. Nevertheless, this trend is much weaker 
than those found by Buat & Xu (1996) and Dale et al. (2007| 
based on the attenuation tracer Lir/Lfuv (cf. Fig. [T3] l instead 
of Apuv- The latter authors have also found an anticorrelation 
of Lir/Lfuv and the specific SFR for intermediate- to late-type 
galaxies, which cannot be confirmed by us. A possible explana- 
tion of these discrepancies could be a systematic weakness of 
Lir/Lfuv. In contrast to Lfuv, ^^ir is related to the obscuration 
of young and old stars. Therefore, significant deviations should 
appear where dust-absorbed light from cool stars significantly 
contributes to Lir (see, e.g., Buat et al. 2005; Cortese et al. 120081 
and references therein). The distribution of the sample galaxies 
in the Afuv-^'Ir/^'Fuv plane shown in Fig.|2T|in comparison to 
the locus of a typical starburst SED for different opacities in- 
dicates that early-type spirals probably exhibit the largest devi- 
ations, while the latest types are close to the starburst curve. In 
detail, the mean deviations of Lir/Lfuv from the reference curve 
amount to 0.31 + 0.07 for SO-Sab, 0.17 + 0.02 for Sb-Scd, and 
0.03 + 0.05 for Sd-I. As indicated by Fig. [21] early-type spi- 
rals also exhibit the largest deviations from the relation between 
i-iR/i-Fuv and Afuv for nearby star-forming galaxies derived by 
Burgarella et al. (l2U05T l (cf. Buat et al. ,2005 ; Cortese et al. 20081. 
Consequently, the apparent trend of Lir/Lfuv with morphology 
and (p could be strengthened/produced by the SFH dependence 
of the luminosity ratio used. 

A reliable trend probably exists for the relation between 
Afuv and Mjtar. The highest obscurations appear to be present 
for masses around 10" M©, i.e. for galaxies in the transition re- 
gion. On the other hand, there is also a relatively large scatter 
in Afuv of high mass galaxies. For low mass galaxies the dust 
obscuration tends to decrease with decreasing mass. Late-type 
galaxies of low mass are obviously characterised by lower dust 
column densities or a higher porosity of the dust screen in front 
of young stars compared to galaxies of earlier type (cf. Dale et 
al. 2007). Apart from irregular and screen-diluting dust distri- 
butions, a "lack" of UV-absorbing dust grains could also be due 
to a deficit of dust production compared to dust destruction or a 
low metallicity. 



4. Discussion and conclusions 

We have developed the SED-fitting code CIGALE as a tool for 
studying basic properties of galaxies in the near and the dis- 
tant Universe. The models constructed by CIGALE consist of 
Maraston ( |2005T l or PEGASE (Fioc & Rocca-Volmerange [T9971 ) 
stellar population models which are reddened by synthetic atten- 
uation curves based on the Calzetti et al. (2000) law and which 
are corrected for spectral lines, and the dust emission templates 
of Dale & Helou (2002). The construction of the SFH by two 
different complex stellar population models with exponentially 
declining SFRs and different amounts of attenuation enables the 
code to deal with age-dependent extinction effects (Silva et al. 
1T9981 Kong et al.[2)04; Panuzzo et al. |2007l Noll et al. |2007] l, 
which is crucial to fit "normal" star-forming galaxies that are not 
characterised by a uniform dust screen in front of the stars (see 
Fig. [T3] l. The models cover the wavelength range from far-UV 
to far-IR, which allows the effect of dust on galaxy SEDs to be 
treated in a consistent way. In the case of dust emission related to 
a non-thermal source, the balance between the stellar luminos- 
ity absorbed by dust and the corresponding emitted luminosity 
in the IR can be preserved by considering an additional hot dust 
component. 
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As our study of the multi-wavelength photometric data of 
a test sample of 39 nearby galaxies selected from SINGS 
(Kenn icutt et al. 12003 1 Dale et al. 120071 Munoz-Mateos et al. 
I2009I I shows, especially the star-formation-related properties 
can be derived with high reliability if the photometry reaches 
from the rest-frame UV to wavelengths greater than lOyum 
(Sect. |3.23T l. Otherwise, the multi -parameter models are not well 
constrained in the IR and the results can be affected by sys- 
tematic errors. In the latter case, the models have to be simpli- 
fied by reasonable a-priori estimates of part of the parameters. 
Apart from the wavelength coverage of the filter set, the individ- 
ual uncertainties in filter fluxes also play an important role for 
the quality of the fit results. Large uncertainties can prevent the 
derivation of parameters that are nearly degenerated regarding 
the shape of the galaxy SED such as details of the SFH or the at- 
tenuation curve. However, it is more critical if too low photomet- 
ric errors overconstrain the model parameters, yielding relatively 
precise but possibly wrong results. Due to unknown systematic 
errors in the object photometry and the models (see Sect. |3.2.2] i 
it is advisable to assume relative errors between 5% and 10% in 
minimum. In contrast, the selected approach for the derivation 
of the expectation values and standard deviations of the different 
parameters does not significantly affect the code results at least 
on average (Sect. I3.2.2I I. The 'max' method introduced by us 
(Sect. |2!2] i and the 'sum' method (e.g., Kauffmann et al. |2003al 
Salim et al. 2007; Walcher et al. I2008I I give similar results as 
long as there are no secondary peaks in the parameter probabil- 
ity distribution which are only poorly populated by models. 

The diagnostics of the SED-fitting results of our SINGS test 
sample has revealed that the most reliable values are obtained 
for non-basic (by the code derived) model properties such as 
the total stellar mass Mstai-, the SFR, the effective age f 04000, 
the bolometric luminosity Lboh the dust luminosity Ldust, and 
the far-UV attenuation factor Apuv (Sect. 13.3. It . Trustworthy 
results are also found for basic input parameters of the mod- 
els such as the burst fraction /burst, the V-band dust attenuation 
of the young stellar population Av ysp, and in part the slope of 
the IR models a. However, the latter properties significantly de- 
pend on the model grid chosen and are not universal, therefore. 
The non-basic properties usually show weak model-related con- 
straints only. An example is the upper fD400o limit of lOGyr in 
the sample due to our restrictions regarding the SFHs investi- 
gated by the code (Sect. [3T2Tt . For the masses given, it has to be 
taken into account that we provide results for Maraston ( 2005 ) 
models and Salpeter IMF (see Sect. l2.1.T] ). In any case, the SFR 
and Ldust derived indicate good agreement with data from other 
studies (Sect lTTTI ). 

An investigation of relations between the different reliable 
galaxy properties has confirmed that the star formation activity 
of nearby galaxies as traced by the specific SFR <p and the ef- 
fective age f D4000 depend on morphology (Sect. 13.3.21 ). Weaker 
trends are also found for Mstar and Apuv- The typical star forma- 
tion activity significantly changes at Mstm- ~ 10" Mq. This mass 
range also indicates the most dust-obscured galaxies in our sam- 
ple. In contrast, far-UV attenuation does not appear to depend on 
star formation activity. These results show what kind of studies 
are possible with CIGALE for the data available for the sample 
of SINGS galaxies investigated. Since the photometric data of 
the SINGS galaxies is characterised by good coverage of the IR 
and comparatively modest coverage and quality in the UV, we 
expect that for other samples especially at higher redshifts the 
reliability of the different model parameters could be different. 
For good quality data in the rest-frame UV and optical, we can 
imagine that stellar population properties and details of the dust 



attenuation could be better studied than it has been possible for 
the SINGS sample. On the other hand, the frequently missing 
information in the far-IR for high-redshift galaxies could signifi- 
cantly affect the quality of Ljust and other star-formation-related 
parameters. However, instruments such as Herschel or ALMA 
will improve the situation in future. Hence, we are convinced 
that CIGALE is a valuable tool for the characterisation of galaxy 
populations in the near and distant Universe. In a series of forth- 
coming papers we will demonstrate this by discussing samples 
of distant galaxies with different selection criteria. 
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